Spin transport in a unitary Fermi gas close to the BCS transition 
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We consider spin transport in a two-component ultracold Fermi gas with attractive interspecies 
interactions close to the BCS pairing transition. In particular, we consider the spin-transport 
relaxation rate and the spin-diffusion constant. Upon approaching the transition, the scattering 
amplitude is enhanced by pairing fluctuations. However, as the system approaches the transition, 
the spectral weight for excitations close to the Fermi level is decreased by the formation of a pseu- 
dogap. To study the consequence of these two competing effects, we determine the spin-transport 
relaxation rate and the spin-diffusion constant using both a Boltzmann approach and a diagram- 
matic approach. The former ignores pseudogap physics and finite lifetime effects. In the latter, 
we incorporate the full pseudogap physics and lifetime effects, but we ignore vertex corrections, so 
that we effectively calculate single-particle relaxation rates instead of transport relaxation rates. We 
find that there is qualitative agreement between these two approaches although the results for the 
transport coefficients differ quantitatively. 



I. INTRODUCTION 

Over the last decade, there has been a growing interest 
in the physics community in the properties of spin trans- 
port, spurred by the idea of using the electron spin as 
a carrier of information. A spin current is a net flow of 
spin and is fundamentally different from a charge current, 
because it relaxes in a different way. In particular, spin 
currents can be strongly affected by interactions, even in 
Galilean-invariant systems. 

Ultracold fermion systems consisting of two spin 
species provide a valuable model system for the study 
of the effects of interactions on spin transport, because 
of their high tunability and the absence of other factors 
which limit the spin conductivity, such as disorder. When 
the cloud of one species moves relatively to the cloud of 
the other species, the latter is dragged along due to mo- 
mentum relaxation by the interatomic interaction, and 
as a consequence, the spin current relaxes. This mech- 
anism is called spin drag pQ. A similar phenomenon, 
Coulomb drag, occurs in bilayer systems, in which the 
drift momentum difference between the carriers in the 
top and bottom layer relaxes due to the Coulomb inter- 
action [SM3]. 

In an important recent experimental study, the spin 
susceptibility, the spin-diffusion constant, and the relax- 
ation rate of spin transport were investigated in a ultra- 
cold two-component fermion gas in the unitarity regime, 
where the interspecies scattering length goes to infinity so 
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that the interactions are as strong as quantum mechan- 
ics allows [H[S]. Inspired by this work, the spin-transport 
relaxation rate and the spin-diffusion constant were cal- 
culated using a Boltzmann approach for high and low 
temperature ranges separately [BJ. 

When an ultracold Fermi gas with two spin species 
and an attractive interaction between the spins is cooled 
to low enough temperatures it shows a transition to a 
superfluid state, where the opposite-spin atoms pair up 
to form Cooper pairs. The effect of interactions on spin 
transport close to this BCS pairing transition is an in- 
teresting subject, since two effects are competing. On 
the one hand, the scattering amplitude between fermions 
is enhanced by pairing fluctuations (preformed Cooper 
pairs) not taken into account in Ref. [BJ. When the tem- 
perature is lowered, this effect ultimately leads to a di- 
verging interaction strength and an instability in the sys- 
tem towards the BCS state at the critical temperature. 
As a consequence, transport coefficients like the spin- 
transport relaxation rate and the spin-diffusion constant, 
are expected to be strongly affected when the system 
approaches the transition. These effects were not seen 
experimentally in Ref. [3], possibly because T c was not 
reached but possibly also due to the competing effect we 
discuss shortly. 

In related work, for ultracold fermions with repulsive 
interactions, enhancement of the spin-transport relax- 
ation rate was predicted close to the ferromagnetic tran- 
sition [7] , and for ultracold bosons close to Bosc-Einstein 
condensation [5]. The latter was recently also observed 
experimentally 9 . Riedl et al. studied the frequencies 
and damping of collective modes in a ultracold Fermi gas 
close to the BCS pairing transition [TU] and found an en- 



2 



hancement of the collision rate in a fermionic gas close to 
the BCS pairing transition. The behavior of the viscos- 
ity was studied in Ref. [TT] . In electron-hole bilayers the 
electrons in one layer and holes in the other can form exci- 
tons, which are expected to condense for low enough tem- 
peratures. Theoretically, it was predicted that the trans- 
port relaxation rate is enhanced when approaching this 
transition 12J. An enhancement was measured experi- 
mentally [IS], although it is still under debate whether 
this enhancement was indeed caused by exciton conden- 
sation. Also for a topological insulator thin film, an en- 
hancement of the transport relaxation rate was predicted 
close to (in this case topological) exciton condensation 

EH- 

However, as already mentioned, there is a competing 
effect at play. For temperatures below the transition tem- 
perature, the system is in the supcrfluid state and the 
excitation spectrum is gapped. Already above T c , a pre- 
cursor of this gap can be seen as the suppression of spec- 
tral weight close to the Fermi level, a so-called pseudogap 
[15] . The suppressed spectral weight is closely linked to 
the reduced lifetime for these excitations. The effect of 
this pseudogap is to reduce the phase space for scattering 
events around the Fermi level, which is expected to re- 
duce the enhancement of transport coefficients. It is the 
competition between the increase in scattering amplitude 
due to pairing fluctuations on the one hand, and the de- 
crease in available phase space due to the pseudogap on 
the other, and its effect on the transport relaxation rate 
and the spin-diffusion constant that we consider in this 
paper. This competition was considered in Ref. [TB] for 
the spin susceptibility. We note related work determining 
the static spin susceptibility for repulsive interactions in 
Ref. [IB]. 

In general, the relaxation rate of a current and the 
decay rate of a single-particle excitation are different be- 
cause of the amount with which the possible scattering 
directions contribute to either case. In the case of the 
decay of a single-particle excitation, all directions are 
weighed equally, but in the case of current relaxation, 
forward scattering (in the direction of the current) con- 
tributes much less than back scattering (in the direction 
opposite to the current). Theoretically, in diagrammatic 
calculations, a transport relaxation rate is obtained by in- 
clusion of so-called vertex corrections to the appropriate 
response function. If these are neglected, the transport 
relaxation time is essentially equal to the single-particle 
relaxation time. 

In this paper we consider a two-component gas of 
fermions consisting of two spin states labeled by a =^, ^. 
We consider the balanced case where the densities of each 
spin component are equal = = n, and where both 
species have the dispersion £(fc) ~ h 2 k 2 /2m — /i. This 
article consists of two parts. In the first part in Sec. [TT] 
we use Boltzmann theory to calculate the relaxation rate 
for spin transport in this system for arbitrary interaction 
strength, both close to the BCS pairing transition and for 
temperatures much higher than the Fermi temperature. 



Then, using the non-interacting spin susceptibility, we 
can determine the spin-diffusion constant using the Ein- 
stein relation. This calculation incorporates the effect of 
pairing fluctuations, but does not take into account the 
pseudogap physics or the effect of the finite lifetime of 
the quasiparticle eigenstates. At the diagrammatic level, 
however, this calculation includes vertex corrections. 

In the second part in Sec. |III| we calculate the fermion 
self-energy at unitarity both close to the BCS pair- 
ing transition and for high temperatures within the 
many-body T-matrix approximation. This self-energy 
is the many-body T-matrix closed with a bare (non- 
interacting) fermion line. Using this self-energy, we deter- 
mine the spectral function, and the spin-transport relax- 
ation rate and spin susceptibility ignoring vertex correc- 
tions. Again, the spin-diffusion constant can be obtained 



using the Einstein relation. We end in Sec. IV with our 
conclusions. 

By comparing the spin-transport relaxation rate and 
diffusion constant obtained using both methods, we as- 
sess the importance of pseudogap physics, finite lifetime 
effects, and vertex corrections. We find that although the 
results of these two methods for the transport coefficients 
differ quantitatively, there is qualitative agreement. This 
indicates that pseudogap physics, finite lifetime effects, 
and vertex corrections do not have a critical influence on 
the prediction using Boltzmann theory for the behavior 
of transport coefficients close to the BCS transition. 



II. BOLTZMANN THEORY 

In this section we derive an expression for the spin- 
transport relaxation rate using the Boltzmann equation. 
We incorporate many-body effects, so-called pair correla- 
tions, close the BCS pairing transition. We apply a spin 
dependent driving force Ff = F and F± = —F which will 
give rise to a spin current j s . The momentum increase 
by the driving force is balanced by the momentum re- 
laxation due to the interaction between opposite spins, 
so that the system reaches a steady state. This mecha- 
nism is called spin drag. In linear response, the quantum 
kinetic equations are written as 



^F-d k n F (ak)) = T t (k) 
lF-d k n F (ak)) = T i (k), 



(1) 
(2) 



where np(e) — 1/[1 -I- exp(/3e)] is the Fermi-Dirac distri- 
bution with j3 = 1/ksT the inverse thermal energy and 
where the collision integral r a (k) gives the net particle 
flux into the (k, \ a)) state due to interspecies interactions. 
The Fermi's golden rule expression for T^(fc) is 



r t (fc) 



2tt 

W 2 



£ki+fe2,fe3+k4 < K£i+£2— (3— £4) 



x |T JM T/T4/i,2(l - A, 3 )(l - hA){5 kl , k - 8 kA . k ), (3) 
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where V is the volume and where we used a shorthand 
notation for the energy £j = £(fej) and for the non- 
equilibrium distributions f a ^ = f a (ki). We note that an 
analogous expression holds for T^(fc). Eq. ^ considers 
the total effect of incoming particles with momenta and 
spin fei,t and k 2 ,i that are scattered into momenta and 
spin fc 3 , J, and k^ 'f . Momentum and energy conservation 
is ensured by the Kronecker and Dirac delta's, respec- 
tively. The many-body T-matrix T MB gives the scatter- 
ing amplitude incorporating the effects of the medium, 
and will be specified below. 

The distribution functions are shifted from the equilib- 
rium distribution by a momentum k a = mv a /h, where 
v is the so-called drift velocity, which is related to the 
spin current by j s = n^v^ — n^v^ — 2nv. Here, we 
used that due to symmetry, the drift velocities of the two 
species will be opposite and of equal magnitude for both 
species. The distribution function f a (k) to first order in 
the drift velocity is f a (k) = «f (£<?(&)) + f£(k), where 
fl{k) = -K{v a ■ k)n' F (^{k)). By subtracting Eq. (B 
from Eq. ([I]), multiplying with hk, and summing over 7c, 
we arrive at the momentum balance equation 



nF=-^ ^(fife)r(fe) =fv + 0(v 2 ), 



(4) 



where the coefficient T of the linear term in the right- 
hand side is given by 



ir(3H 
6^3 



<5fei+fe 2 ,fe3+fc4 ( 5(£i + £2 — & - £4) 



x \T MB \ z nm 2 (l - n 3 )(l - n 4 )(fei - k 2 + k 3 - k 4 y 



(5) 



and we introduced the shorthand notation rij = 
np(£(fcj)). We introduce the spin drag conductivity 
od via j s = <7£>F and from Eq. Q we identify op, = 
— 2n 2 /r. The spin drag conductivity can be related to 
the spin-transport relaxation rate 1 /Yd by the Drude for- 
mula op, = 2nTu/m leading to 1/rp, = —T/mn. 

To make further progress, we rewrite Eq. |5]) in a form 
that is convenient to account for a divergence in the pair- 
ing channel due to pairing fluctuations close to the super- 
fluid transition. First, we introduce an additional energy 
integral over the variable tko — £(fci) + £(^2), i- c -, the 
sum of the energies of the two incoming particles in the 
scattering event. This choice is convenient because close 
to the superfluid transition, in the on-shell approxima- 
tion, the dominant energy on which T MB depends is hui. 
We stress that this choice is different from the conven- 
tional one, where hui is taken to be the difference between 
the energies of the incoming and outgoing ^-particles [3] . 
Next, we resolve the momentum conserving Kronecker 
delta £fe 1+ k 2 ,fe 3+ fc 4 in Eq. (5| by choosing k x = k + K/2, 
k 2 = K/2 -k,k 3 = k' + K/2, and fc 4 = K/2 - k' and 



introduce an additional energy integral via 

<5(£i+£ 2 -£3-£ 4 ) = J d(M<5(£i+£ 2 -M^(£3+£4-M- 

(6) 

Then, after using the Dirac identity Qf[l/(a: + iO)] = 
—ttS(x), we rewrite Eq. ^ into 



1 

T~D 



TT/3H 



dk dk' dK 



6ir 2 mn J (2tt) 3 (2tt) 3 (2tt) 3 
x \T MB (K,hu)\ 2 
sinh 2 (/3/iu;/2) 



d{huj){k' + k) 



x is 



1 - n F (£(fe + K/2)) - n F {Z(K/2 ~ k)) 



huj + »0 - (£(fe + K/2) + ((K/2 - k)) 
1 - n F (g(fc' + K/2)) - n F (((K/2 - fcQ) 
Kuj + »0 - (£(fe' + K/2) + ({K/2 - k')) 

(7) 

The effective interaction is given by the many-body ma- 
trix element T MB (K, hui) which takes a simple form close 
to the superfluid transition in these variables 



T MB (K,huj) 



V 



l-V E(K,huj)- 



(8) 



where Vq = ^itah 2 /m with a the s-wave scattering length. 
The bare pairing susceptibility 5 is given by 



k 

l-n F (£(fe + K/2)) 



n F (£(fe - K/2)) , f 



huj + iO-(((k + K/2)+((k-K/2)) 2e(k) 



(9) 



with e(k) — h 2 k 2 /2m. The second term of the sum- 
mand cures the ultraviolet divergence of summing the 
first term. Its derivation can be found in Ref. [20]. We 
note that the transition temperature T c to the superfluid 
state is determined by I/T MS (O,0) = 0. In Eq. Q, we 
can perform the integrals over k, k' , and the angles of K 
exactly. The integrals over u> and the length of K must 
then be performed numerically. Our theory does not take 
into account the superfluid gap or the presence of Cooper 
pairs for T < T c , and is thus only valid for temperature 
higher than the transition temperature. 

To obtain results for 1/td at constant density, we 
need to solve the non-interacting equation of state n = 
(1/V) Y^,k n iK£(k))i which yields the function /j,(T,n). 
For low temperatures, /i is given by the familiar relation 
fx = €f[1 — (7r 2 /12)(fc_BT/e_F) 2 ]. Here, the Fermi energy 
ep = kpTp = h 2 k 2 F /2m with the Fermi wavenumber 
kp = (67T 2 ?!) 1 / 3 . The transition temperature can then be 
determined as a function of the dimensionless interaction 
parameter kpa. We find at unitarity (—l/kpa — 0) T c 
0.50Tf and that for weak coupling where (—1/kpa 1) 
T c (xTpexp{-n/2kp\a\). In Fig. [I] we show the dimen- 
sionless spin-transport relaxation rate h/epTo obtained 
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FIG. 1: The dimensionless spin-transport relaxation rate 
H/efTd obtained from Eq. |7| as a function of the reduced 
temperature T/Tf- From top to bottom the curves corre- 
spond to interaction parameters -l/fc^a = {0, 0.56, 1.0, 1.4}. 
Each curve terminates at the value of the transition temper- 
ature corresponding to its fei value: in decreasing order 
these are T/T F « {0.50,0.24,0.13,0.064}. The dashed line 
is the high-temperature asymptote 1/td oc l/VT predicted 
in Ref. [B]. 



from Eq. ([7]) as a function of the reduced temperature 
T/Tp. Each curve terminates at the value of the transi- 
tion temperature corresponding to its fcpa value. These 
curves illustrate a number of effects. Firstly, we see that, 
for all fcj70, the relaxation rate is enhanced (but remains 
finite) when approaching T c from above. This enhance- 
ment is larger in the strong-coupling case than in the 
weak-coupling case. The reason for this enhancement is 
that pairing fluctuations become important close to the 
superfluid transition, which lead to the appearance of a 
pole in T MB in the integrand of Eq. Q at K = and 
u = 0. 

For large temperatures, all curves approach 
the dashed line in Fig. [T] given by 1/td = 
(e F /^)(32 % /2/97r 3 / 2 ) v /T F /T which falls off as l/VT, 
and was determined in [5]. The 1/vT dependence can 
be understood by noting that 1/r oc nvcr, where n is 
the density, v the mean particle velocity, and a the 
cross section at unitarity. For temperatures T ^> Tp, 
the scattering cross-section is given by the square of the 
thermal de Broglie wavelength, and thus a decreases like 
1/T. Since for high temperatures it holds that v oc VT, 
we find 1/r oc 1/y/T. This high-temperature behavior 
was also measured experimentally [3]. 

In the weak-coupling case (—1/kpa > 1), we note an 
increase of 1/td with temperature up to T w Tp. An in- 
crease in temperature leads to a larger phase space avail- 
able for scattering events due to reduced Pauli block- 
ing. In the absence of many-body effects (T MB = Vq), 
this effect would lead to the standard T 2 dependence of 
transport coefficients in a Fermi liquid. Indeed, when 
rpMB _ we recover the weak-coupling result by 
Bruun, 1/t d = {16ir/9)(k 2 F a 2 e F /h){T/T F ) 2 for T«T C 
i- 



We note that our theory overestimates the T c at uni- 
tarity in comparison with the value T c = 0.15T F obtained 
by Monte Carlo and RG methods (see e.g. Ref. [TH] ). In 
this sense, the experimental result for the spin-transport 
relaxation rate in Ref. [1] should not be compared to 
the top curve in Fig. [T] at unitarity, but instead to the 
second from below, which has a finite kpa value and 
T c = 0.13Tf- Then, we see that our result agrees quali- 
tatively with the experimental results in Ref. 0]. 

To obtain the diffusion constant D s we use the Einstein 
relation D s = (Td/Xs- For the spin susceptibility, we 
use the non-interacting result Xs = —(1/V) J2k n '(£(k)). 
The result is shown in Fig. [2] where we plot the diffu- 
sion constant scaled by m/H versus the reduced temper- 
ature. The vertical order of the curves is reversed with 
respect to Fig. [I] Each curve again terminates at the 
value of the transition temperature corresponding to its 
kpa value. The remarks made about Fig. [TJfor a large 
part carry over to our results for the diffusion constant 
in Fig. [2] We see a suppression of D s when approach- 
ing the transition for all kpa values. For the top two 
curves with — 1/kpa = {1.0,1.4} we see that the non- 
monotonic behavior of 1 /td is also present in the curves 
for D Sl which for these k F a values now have a mini- 
mum around T/T F = 0.5. Since Xs follows the Curie 
law Xs = n/ksT for high temperatures , we find that for 
high temperatures D s = (h/m)^' 2 /lQy/2){T /Tpf 2 , 
which is shown as the dashed line in Fig. [2] We note the 
factor 2 difference with the result for this asymptote from 
Ref. [6], which is caused by a different definition of op>- 
We note that if we again compare the second curve from 
above, which has a finite kpa value and T c = 0.13Tfj 
with the experimental results from Ref. [4j at unitar- 
ity, we see qualitative agreement: both curves flatten off 
for temperatures below T /Tp — 0.5. Considering our re- 
sults, we see that this behavior is actually the crossover 
behavior between the monotonic behavior of the unitary 
(bottom) curve in Fig. [2] and the non-monotonic behav- 
ior of the weak-coupling (top) curve. Clearly, a strength 
of our approach is that we can determine the whole tem- 
perature range using a single approach. 

III. DIAGRAMMATIC THEORY 

In the second part of this work, we use a diagrammatic 
approach to determine the self-energy, spectral functions, 
and transport coefficients. 

A. Definitions 

First, we recall some aspects of the microscopic theory 
of an interacting gas of fermions. The central quantity in 
this theory is the self-energy, which is the sum of all one- 
particle reducible diagrams, and the spectral function de- 
rived from it. For reference we give the expressions for 
the spectral function, and the density of states and the 
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FIG. 2: The dimensionless spin-diffusion constant D s m/h 
obtained from the result of Fig. [I] using the Einstein relation. 
From bottom to top the curves correspond to interaction pa- 
rameters — 1/kpa — {0, 0.56, 1.0, 1.4}. Each curve terminates 
at the value of the transition temperature corresponding to its 
k F a value. The dashed line is the high-temperature asymp- 
tote D s = {h/m)(9-K 3/2 /l6V2)(T/T F f /2 
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FIG. 3: The self-energy E is the many-body T-matrix closed 
with a single non-interacting Green's function line. The ex- 
ternal lines are indicated by the dashed lines 



particle density in an interacting system in terms of the 
spectral function. For a system which is translationally 
invariant in space and time, the retarded fermion Green's 
function is given by 

-h 



G(k,w) = 



(10) 



where E(fc,w) is the retarded self-energy. The spectral 
function describes the weight of a single-particle excita- 
tion and is given by 



p(k, ui) 
1 



1 



3[G(fc,w)] 



nh [-w + £(fc)//i + 3?E(fc,w)] 2 + 3 2 E(fc,w)' 



(11) 



In the absence of interactions, E = and the spec- 
tral function is a <5-peak at Huj — £(fc). A nonzero 
real part of E will shift the 5-peak to the solution of 
uj = ^(k)/h + !RE(fc,u;), while a nonzero 3E will lead 
to a broadening of the peak, indicating a finite lifetime 
of the single-particle excitation. In particular, a con- 
stant 3E = — 1/2t will lead to a spectral function with a 
Lorentzian shape with width r. We note several relations 
for the spectral function that we use later. The sum rule 
states that 



J dhujp(k,u}) = 1, 



(12) 



and follows from the anticommutation relation of the 
fermion fields. The density of states is 

1 



(2tt) 3 

and the equation of state is 



dkp(k,Lo), 



n = / dhujnF{uj)v{uj). 



(13) 



(14) 



Using the Kubo formalism, we obtain the spin-transport 
relaxation rate 1 /td and spin susceptibility Xs = d(n^ — 
n^)/d(fj,-f — p\J) in terms of the spectral function. These 
are in the absence of vertex corrections given by 



1 



and 



3mn 



— fe 2 / dhuip 2 (fc,w)ri' (huj) 

7- •> 



(15) 



V 



dhcudhuj' 



x p(k,u)p(k,u') 



np(huj) — np(huj') 
hu — tvjj' 



(16) 



We note that Eq. |15|) can be derived by calculating the 
ordinary Drude conductivity of the atoms with spin up, 
say, in the presence of scattering from an 'external' poten- 
tial which is due to the presence of atoms with spin down. 
We note that for E = —i/2r, Eq. (TTsl) gives 1/t d = 1/r, 



and that in the non-interacting case for low tempera- 
tures, Eq. (16 1 reduces to the well known \o — 3n/2eF- 
Finally, we note that we can obtain the spin-diffusion 
constant D s via the Einstein relation D s = cfd/Xs- 



B. Self-energy 

We calculate the fermion self-energy which is appropri- 
ate for this system both close to the BCS pairing transi- 
tion and for high temperatures at unitarity [31]. Indeed 
the many-body T-matrix closed with a non-interacting 
fermion line shown in Fig. [3] incorporates the effect of 
enhancement of the interaction close to the BCS tran- 
sition. In the dilute limit for weak interactions, where 



G 



rpMB _^ ft rec [ uces ^ Q £ ne standard mean- field shift 
of the dispersion VqTi. The expression for £ as a function 
of momentum k and fermionic Matsubara frequency u n 
is 



S(fe, ifuj n ) = 



—- y 



x T MM {K, ihw m )G a {K - k, ih(u m - un)), (17) 

where the summation is over the bosonic Matsubara fre- 
quencies uj m and Go is the non-interacting Green's func- 
tion. We compute the Matsubara summation of Eq. ( 17) 
to obtain 



1 1 

nv 



J2[nF(Z(K-k))+n b (t(K-k)+hw)\ 



K 



x^s[T MB {K,i{K -k)+Hu)]. (18) 

Using $[T MB (K,Hu/)\ = ^[~{K,hu)]\T MB {K,huj)\ 2 , 
we can rewrite Eq. (fl8| as 



3?£(fc,/iw) 



7T 

W 2 



E 



x 6(hw + ^(k 2 ] 



C( fc 3) - £(fc4))5fe+fe 2 ,fe 3 +fe4 



x \T MB (k + k 2 ,hu + t(k2W 
x [n 2 (l - n 3 )(l - rii) + (1 - 712)713714], (19) 

where we again used the shorthand notation m = 
ni?(£(fcj))- Since is a single-particle relaxation rate, 
we see that this rate is exactly the sum of the in and out 
rates for the state with momentum k and energy ftw as 
would be obtained from a Fermi's golden rule expression. 
In Eq. (19), we introduce the center-of-mass momentum 



K = k + fe 2 , and an additional energy integral 

*(fod + £(fc2)-£(fcs)-£(fe4)) 

= J d{W)8{tko + £(fc 3 ) - hn)5(t(k 3 ) + c(fe 4 ) - h£l). 

(20) 

After some rewriting, we then arrive at the following ex- 
pression for 



n^(w)(l — 7lF(w))SJ£(fe, Huj) 



8ir 2 kh 3 



d(hn)dK 



R 



\T"»{K,m)\* 
sinh 2 (/3fifi/2) ' n 

x [1 — n F (hVL — kuj) — n F (huj)}, (21) 

which is a convenient starting point for numerical eval- 
uation. In going from equation Eq. (19) to Eq. (21), we 
introduced OS from E q. pi ). The delta functions on the 
right hand side of Eq. (|20[) restrict the integration range 
to R which is defined by the following inequalities 

2 f i + huj- -e(K) > and fx + W. - Huj > (22) 
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FIG. 4: The imaginary (solid line) and real (dashed line) 
parts of the self-energy at unitarity for hk — -^/2m/i for 
k B T/n = 1 (a) and T = T c (b). 



and 



\y/e(Je)- y/fi + hn-hu}\ < y/e(K) 

< Ve(fc) + V/i + hn-hiu, (23) 

where we defined e(fc) = h 2 k 2 /2m. Then, the real part 
of the self-energy is obtained by a Kramers-Kronig trans- 
form 



,3?S(fe,Sw) 



(24) 



In the next subsection, we show the results for the spec- 
tral function and spin-transport coefficients obtained by 



evaluating Eqs. (21) and (24) 



Results 



By evaluating Eqs. (21) and (24) numerically, we obtain 
Fig. [4] where we show the imaginary (solid line) and real 
(dashed line) parts of the self-energy for a typical momen- 
tum hk — y/2mfi at unitarity. Part (a) of Fig. [4] shows 
the results far from the transition at T = \xjkB ~ 1-5T C , 



7 



pp 





0.30 








A 0.25 




\T=T 

I r 




/ /Y\ 0.20 








// \ 0,15 








/ V- 10 " 








Jr = u/k B oi)5 











fioj/p 



FIG. 5: The spectral function p at unitarity for hk = \/2m/i 
for fcaT/^i = 1 (dashed line) and T = T c (solid line). 



and part (b) of Fig. 
T — T c . From Fig. 



4 1 shows the results at the transition 



we see that upon approaching the 
transition, a divergence develops at uj = in the imagi- 
nary part of the self-energy and a corresponding discon- 
tinuity in the real part. Closer inspection shows that the 
imaginary part diverges logarithmically. When di- 
verges as x log(Aw), with Aw the distance from the pole, 
the resulting jump in the real part is nx. For general fc, 
the position of the divergence is hu = — £(fe), which can 
be understood as follows: the transition is accompanied 
by a pole in the many-body T-matrix T MB (K : Ml) for 
K = and Ml — 0. In the expression for Eq. ( 19 ) 
the T-matrix enters as T MB (k + k2, hio + so that 

for fno = — £(fe), the conditions K = and Ml = are 
satisfied when k 2 = —k. 

As a check, we may also obtain the logarithmic di- 
vergence analytically. The divergence is caused by the 
development of a pole at Vl = and K = in Eq. (21 1. 
For small fi and K, T MB cx l/{a 1 n+a 2 K 2 +a 3 (T-T)j+ 
ia^Vi) (with the a,'s real) while the rest of the integrand 
is linear in K. Integration over f2 and K then yields a 
divergence log(T — T c ). In contrast, for small fl and K 
the integrand in the Bolzmann case Eq. Q goes as K 2 . 
Integration shows that there is then no divergence, as 
was shown in Fig. [T] 

From Eq. ( [IT] ) we see that a divergence of will lead 
to a suppression of the spectral function p. Indeed in 
Fig. [5] we show the spectral function p at unitarity for 
hk = y / 2to/z, for ksT/p = 1 (dashed line) and T = T c 
(solid line) . Already quite far from the transition at T = 
pt/ks ~ 1-5T C , interactions cause a suppression of the 
spectral weight close to u) — 0, which moves into peaks 
on both sides of lo — in order to keep the sum rule 
satisfied. At T = T c , this suppression is maximal, and 
here p vanishes logarithmically as \lu\ —> 0. 



In Fig. [6] we show the solution of the equation of state 
n = J d{hui)nF{ui)v{io) for \i. The dashed line corre- 




-1.5 L 



FIG. 6: The solution of the number equation giving p as a 
function of temperature. The dashed line is the non interact- 
ing result. The solid line gives the result at unitarity. 
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FIG. 7: The density of states v versus u>. The dashed curve 
is the non-interacting result for T = 0, while the solid curve 
correspond to the unitary case for T — T c . 



sponds to the non- interacting result, which was used to 
obtain the Boltzmann results in Fig. [JJ The solid line is 
the solution of the number equation at unitarity, which 
terminates at the critical temperature T c w 0.29Tp. The 
difference of this value for T c with the mean-field value 
for T c at unitarity T c ps 0.5TV used in Sec.|TTJis that there, 
we used the non-interacting equation of state. That 
p is lower than the non-interacting value is intuitively 
clear: an attractive interaction will lower the energy of 
the atoms and thus p must also be lowered in order to 
keep the number of particles the same. 

In Fig. [7] we show the density of states v versus uj. The 
dashed curve is the non-interacting result at T = 0, which 
goes as v oc -s/uj — ep. The solid curve is the result at 
unitarity at T = T c . We see that the suppressed spectral 
weight around to = leads to a greatly reduced density 
of states, a so-called pseudogap. We note that since T c 
is not particularly low, the spectral functions are never 
very sharply peaked, and the density of states remains 
nonzero in the pseudogap in agreement with the results 
of Ref. [T5]. 
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FIG. 8: The spin-transport relaxation rate 1/td as a func- 
tion of temperature at unitarity. The dashed line is the Boltz- 
mann result, the solid line the diagrammatic result obtained 
by evaluating Eq. ( 15 1. 
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FIG. 9: The spin-diffusion constant D a as a function of tem- 
perature at unitarity. The dashed line is the Boltzmann re- 
sult, the solid line the diagrammatic result obtained by evalu- 
ating the spin susceptibility in Eq. ( 16 1 and using the Einstein 
relation. 



In Fig. |8] we show the spin-transport relaxation rate ob- 
tained by evaluating Eq. ( 15 ) at unitarity as the solid 
line. The dashed line is the unitarity result from the 
Boltzmann calculation (the top line in Fig. [IJ). Since 
both lines are determined using different equations of 
state, they terminate at different T c 's. Similarly, in Fig. [9] 
we show the spin-diffusion constant. The dashed line is 
the unitarity result from the Boltzmann calculation (the 
bottom line in Fig [2]), while the solid line is obtained by 
evaluating the spin susceptibility in Eq. ( 16 ) and using 
the Einstein relation to find D*. 



In Fig. [To] we show the spin susceptibility \s obtained 
by evaluating Eq. ( 16 ) scaled by the non-interacting zero- 



temperature value xo = 3n/2e^ as the solid line. The 
dashed line corresponds to the non-interacting result. We 
see that for large temperatures the lines converge, and 
that our result shows a downturn close to T c . This down- 
turn is expected physically, since the magnetic response 
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FIG. 10: We show the spin susceptibility obtained by evalu- 
ating Eq. ( 16 1 scaled by the non-interacting zero-temperature 
value xo = 3n/2e F as the solid line. The dashed line corre- 
sponds to the non-interacting result. The solid curve termi- 
nates at the critical temperature T c = 0.29Tf. 



should diminish, when f and J, spins become more corre- 
lated [TB]. We note that in BCS theory the downturn of 
the spin susceptibility happens below T c . 



IV. CONCLUSIONS AND DISCUSSION 

We have determined the spin-diffusion constant and 
the spin-transport relaxation rate using a Boltzmann ap- 
proach and a diagrammatic approach at unitarity, as 
shown in Figs. [8] and [9] On a diagrammatic level, the 
Boltzmann calculation takes into account vertex correc- 
tions, while our diagrammatic approach does not. Op- 
positely, the diagrammatic approach takes into account 
pseudogap physics which suppresses the spectral weight 
close to the Fermi level, while the Boltzmann approach 
does not. Surprisingly, we find qualitatively equivalent 
behavior for both transport coefficients using these ap- 
proaches. Seemingly, pseudogap physics and vertex cor- 
rections are not of critical importance when evaluating 
these transport coefficients close to the BCS transition. 
To research this claim further, vertex corrections appro- 
priate to the self-energy we used should be evaluated. 

We note a recent preprint containing related work 22J , 
in which Luttinger-Ward theory is used to study spin 
diffusion in Fermi gases. 
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